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1 0 Related Applications 

This patent application is a continuation-in-part of U.S. patent application serial 
no. 09/575,492, filed May 18, 2000, entitled Method and Apparatus for Determining 
Global Position Using Almanac Information, which is hereby incorporated by reference 
in its entirety. 

1 5 Field of Invention 

The invention relates to satellite based positioning methods and systems, and 
more particularly, to using stored Global Positioning System (GPS) satellite almanac 
information to determine the coarse position of a mobile station. 
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Background 

Satellite based positioning systems include networks of earth orbiting satellites 
. that constantly transmit geographical position information to receivers. An example of a 
satellite based positioning system is the Global Positioning System (GPS), which 
5 includes a network of earth orbiting satellites (GPS satellites or satellite vehicles). A 
GPS receiver can be in a fixed position on a base station or location server, or the 
receiver can be on a mobile unit or station. The GPS receiver and known processing 
hardware and software are used to receive GPS data from GPS satellites and determine 
the position of the receiver. GPS technology was first used in military and emergency 
10 service applications. As GPS technology becomes more economical and compact, 
however, it is becoming ever more common in consumer applications. For example, 
many location based services (LBS) now exist, such as asset tracking, turn-by-turn 
routing, and friend finding. 

The wireless community is investigating procedures and methodologies for LBS 
1 5 that most benefit both the consumer and the provider. The change in focus from 

emergency to commercial services is characterized by an increased number of users, an 
increased number of positioning requests per user, and an increase in the accuracy needed 
to meet various LBS quality of service requirements. As the demand for LBS capabilities 
grows, the number of users taking advantage of those services will potentially outstrip 
20 those transmitting requests for emergency services. Additionally, a number of these 

services (turn-by-turn routing, for example) involve each user requiring repeated location 
solutions, which is in contrast to the single position request that characterizes an 
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emergency call Finally, many of the LBS that are being planned require accuracies far 
greater than those required for regulatory compliance. 

The diversity inherent in proposed LBS offerings presents LBS providers with 
many choices and challenges. For example, some LBS work better with a mobile station 
5 (MS) based approach while others are better served by an MS-assisted solution. 

It is desirable to have a satellite based positioning method and system that is able 
to handle greatly increased LBS system usage while allowing a system operator the 
flexibility to choose between MS based, MS-assisted, or other approaches. It is further 
desirable to have a satellite based positioning method and system with the high accuracy 
1 0 required by most LBS. It is desirable to have a satellite based positioning method and 
system that minimizes the use of network resources while also minimizing the power 
used by mobile stations. 

Brief Description of the Drawings 

The following figures are provided to assist in describing embodiments of the 
15 invention, and are not intended to be exclusive or limiting. In the figures, like reference 
numerals designate like elements. 

Figure 1 is a block diagram of an embodiment of a satellite based positioning 

system; 

Figure 2 is a flow diagram of a general embodiment of a coarse location 
20 positioning method; 
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Figure 3 is a flow diagram of another embodiment of a coarse location 
positioning method; 

Figure 4 is a flow diagram of another embodiment of a coarse location 
positioning method; 

5 Figure 5 is a flow diagram of another embodiment of a coarse location 

positioning method; 

Figures 6A through 6H illustrate the results of one simulation of an 
embodiment; and 

Figures 7 A through 7H illustrate the results of another simulation of an 
10 embodiment. 
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Detailed Description 

Embodiments of the invention allow system operators providing LBS to flexibly 
perform position fixes using various combinations of network resources and mobile 
station resources. Embodiments of the invention allow system operators to provide LBS 
5 to large numbers of subscribers with high accuracy within acceptable time frames. 

Embodiments of the invention allow coarse positioning by a mobile station in a 
satellite based positioning systems, such as the GPS system. Embodiments of the 
invention provide traditional GPS accuracies while dramatically reducing the bandwidth 
consumed by similar methods and their required assistance data payloads. Embodiments 

10 of coarse location positioning described herein require a minimum amount of transmitted 
assistance data. While some MS-assisted techniques rely on acquisition aiding 
information with a usable lifespan measured in minutes, or on ephemeris data with a 
lifespan measured in hours, coarse location positioning methods described herein use 
satellite sub-almanacs that remain usable for many weeks. An almanac is a list of Kepler 

15 parameters, e.g., orbit and clock parameters, for all of the satellites operating in the GPS 
constellation. A sub-almanac contains only the orbit and clock parameters for one 
satellite. 

Complete GPS almanac information is assembled from data broadcast by each 
satellite in a twelve-second message every thirty seconds, for a total time of twelve 
20 minutes. A sub-almanac can be retrieved by receiving any twelve-second message within 
the twelve minute broadcast. Each sub-almanac is identified by a unique sub-almanac 
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identification. The Kepler parameters are typically not considered accurate enough to be 
used for precise positioning calculations. 

Each satellite continuously transmits its own ephemeris data. Ephemeris data 
contains satellite position information that is accurate to within five to ten meters of the 
5 actual satellite actual location. However, the ephemeris data is only valid for about four 
hours. The ephemeris data is broadcast from each satellite for eighteen seconds each 
thirty-second period. 

Embodiments described herein allow mobile stations to perform "off-line" 
processing in which a mobile GPS receiver that has no communication with a network, 

10 even for several days, can calculate its own coarse position using stored almanac data. 
This allows the use of significantly less power than conventional methods. The mobile 
station can communicate the coarse position, and identifications of sub-almanacs used in 
the coarse position calculation, to a server on the network when communication is 
restored. At that time, corrections to the position can be made with the help of the 

15 network server. In various embodiments, corrections to the position are made at the 
pseudorange level or at the three-dimensional position vector level with a three- 
dimensional position correction vector. For example, the corrections can be calculated as 
a pseudorange correction per satellite. In this case, the sub-almanac based pseudoranges 
are extracted back from the coarse position, and the pseudorange corrections are applied 

20 before recomputing a position with the ephemeris and with accurate time. In another 
example, the corrections are calculated as a position vector that is calculated over the 
same satellites used to compute the coarse position. 
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Figure 1 is a block diagram of a GPS system 100. The system 100 includes a 
mobile station 140 (MS) that communicates with a network 1 10. The network 1 10, in 
one embodiment, includes a portable or fixed base station (BS) 130, and a portable or 
fixed location server 120. The location server 120 acts as a data server to the BS 130, 
5 which in turn communicates with the MS 140. The location server 120 has access to 
GPS data from satellites in view. Typically the location server 120 stores the current 
ephemeris and almanac data for at least the satellites in view. In the embodiment shown, 
the location server 120 includes a GPS receiver 121 for receiving the GPS data directly. 
In other embodiments, the location server does not include a GPS receiver, but receives 
10 the GPS data transmitted from another component that has a GPS receiver. The location 
server 120 further includes a central processing unit (CPU) 122, a memory device 123, 
and a communications channel interface 124. The communications channel can include 
any type of communications channel, such as wired or wireless. 

BS 130 includes a central processing unit (CPU) 131, a memory device 132, a 
15 wireless transceiver 133, and a communications channel interface 134. 

The MS 140 includes a GPS receiver 141, a CPU 142, a memory device 143 and a 
wireless transceiver 44. The MS 140 has sub-almanac data stored in the memory 143, 
although the data may not be a complete almanac. 

In other embodiments, the BS 130 also includes GPS capability. Alternatively, 
20 complete GPS capability and all of the capability necessary to communicate with the MS 
140 is collocated. In general, the MS 140 communicates with some configuration of 
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components on the network 110. The network 1 10 provides LBS such as asset tracking, 
in which mobile assets (e.g., taxi cabs) are mobile stations with GPS receivers. 

Embodiments of the invention reduce the time to first fix (TTFF) for a GPS 
mobile station receiver while using minimal power, in part by using GPS sub-almanac 
5 data stored on, or sent to, the mobile station to calculate its own coarse position. Figure 
2 is a flow diagram illustrating a general method according to one embodiment. At 200, 
a coarse positioning session is initiated. The coarse positioning session can be initiated 
by the MS itself, or the MS may receive a signal initiating the session, such as a 
positioning request sent or forwarded by the position server or BS. Optionally, the 
10 network also transmits a reference position that includes latitude, longitude, and time. 
The reference position is a precise position of a location that is close enough to the 
mobile station to be helpful to the mobile station in as an indication of which satellites 
should be in view and in finding an initial position solution. 

At 202, the MS uses stored sub-almanac data along with a reference position (if 
15 received) and time estimates, to acquire satellites and take measurements. Taking 

measurements includes, for example, determining a approximate distance of the mobile 
station from a particular satellite with some error attributable to time differences. The 
MS, at 204, uses the measurements, and a satellite position as derived from the sub- 
almanacs, to calculate a coarse position. At 206, the MS transmits the coarse position to 
20 the network along with a list that identifies the particular satellites used in the solution, 
and the particular sub-almanacs used for each satellite. In some embodiments, the MS 
transmits a full coarse position. In other embodiments, the MS transmits a position 
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difference between the reference position received from the network, and the coarse 
position. 

At 208, the network derives satellite positions using the current ephemeris. The 
network then calculates the differences between almanac-derived and ephemeris-derived 
5 satellite positions and uses them to create a corrected position for the MS at 210. The 
network, in various embodiments, can calculate a coarse correction in various ways. For 
example, the correction can be a pseudo-range correction per satellite. In this case, the 
sub-almanac pseudoranges are extracted back from the coarse position, and the 
pseudorange corrections are applied before recomputing the position with the ephemeris, 

10 and with time. Alternatively, the coarse position correction can be calculated as a 

position correction vector that is calculated over the same satellites as those used by the 
MS to compute the coarse position. The coarse position correction may or may not 
include differential corrections in different embodiments, where differential corrections 
account for discrepancies between what the ephemeris predicts and what a reference 

1 5 receiver at a known location sees in the pseudoranges. 

At 212, the network determines whether one or more of the sub-almanacs used by 
the MS contributed an unacceptable level of range error, due for example to the sub- 
almanacs being too old. If so, the network transmits the latest version of the particular 
error-causing sub-almanacs to the MS at 214. 

20 A sub-almanac stored at the MS is only replaced when the error associated with 

the aging data exceeds a specified threshold. This corresponds to a useable life for a sub- 
almanac that is approximately three weeks to twenty weeks. In one embodiment, an 
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indication of the acceptable level of error is transmitted from the MS to the network. The 
MS uses the sub-almanacs for both satellite acquisition and for calculating a coarse 
position. 

In various embodiments, the order of the numbered process elements shown in all 
5 of the figures, including Figure 1, may be changed. Alternatively, all of the numbered 
process elements shown may not be present in various embodiments. For example, in 
some embodiments, the mobile station may calculate and store a coarse position using 
stored sub-almanacs, even though the mobile station may not be in contact with the 
network for some period of time, such as several days. When network contact is restored, 
1 0 the mobile station can immediately transmit its coarse position. 

The coarse location positioning method illustrated in Figure 2 may be applied 
using a variety of specific schemes. Figure 3 is a flow diagram illustrating one 
embodiment of the coarse location positioning method. 

At 300, a coarse positioning session is initiated. Optionally, the network also 
1 5 transmits a reference location that includes a latitude, longitude, and time. At 302, the 
MS uses stored almanac data along with its local position and time estimates, to acquire 
satellites and take measurements. The MS, at 304, uses the measurements, and satellite 
positions as derived from the sub-almanacs, to calculate a coarse position. At 306, the 
MS transmits the coarse position to the network along with a list that identifies the 
20 particular satellites used in the solution, and the particular sub-almanac used for each 
satellite. 
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At 308, the network calculates an estimated range error for each sub-almanac. At 
3 1 0, the network determines whether any of the sub-almanacs used by the MS exceed the 
range error threshold. If none of the sub-almanacs used by the MS exceed the range error 
threshold, the network calculates the final position solution at 312. If one or more of the 
5 sub-almanacs used by the MS exceed the range error threshold, the network transmits 
replacement sub-almanacs for each of the sub-almanacs exceeding the range error 
threshold at 3 14. The network then reissues the position request at 3 16. In an LBS 
service characterized by repeated position requests, the single- request path (the path 
including block 3 12) is traversed in the vast majority of sessions. 

10 Figure 4 is a flow diagram illustrating another embodiment of the coarse location 

positioning method. At 400, a coarse positioning session is initiated. Optionally, the 
network also transmits a reference location that includes a latitude, longitude, and time. 
Upon receipt of the position request, at 402, the MS transmits a list of the satellites it will 
track and the associated sub-almanac identifications to the network. At 404, the network 

1 5 calculates the estimated range errors for each sub-almanac and transmits replacement 
sub-almanacs for those exceeding the range error threshold. At 406, the MS uses its 
stored sub-almanacs, including any received replacement sub-almanacs, to acquire 
satellites and take measurements. At 408, the MS calculates a coarse position using the 
sub-almanacs. At 410, the MS transmits the coarse position to the network along with a 

20 list of the satellites used in the solution and the associated sub-almanac identifications. 

The network calculates the final position solution at 412 from the information transmitted 
from the MS. 
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Figure 5 is a flow diagram illustrating another embodiment of the coarse location 
positioning method. At 500, a coarse positioning session is initiated. At 502, the MS, 
upon receipt of a position request, calculates which satellites it will track. At 504, the 
MS examines the age of the sub-almanac data associated with each satellite in its 
5 proposed tracking list and determines whether the age of any of the sub-almanacs is 
above a value derived from the predetermined range error threshold. If none of the sub- 
almanac ages are above the value derived from the range error threshold, as shown at 
506, the MS uses its stored sub-almanacs to acquire satellites and take measurements. At 
508, the MS calculates a coarse position using the sub-almanacs. At 510, the MS 
10 transmits the coarse position to the network along with a list of the satellites used in the 
solution and the associated sub-almanac identifications. The network calculates the final 
position solution at 512 from the information transmitted from the MS. 

If any sub-almanacs exceed an age threshold, as shown at 514, the MS transmits 
the list of the satellites it intends to track, the satellites' associated sub-almanac 

1 5 identifications, and the range error threshold to the network. At 5 1 6, the network 

calculates the estimated range errors for each sub-almanac and transmits replacement 
sub-almanacs for those exceeding the range error threshold. At 5 1 8, the MS uses the 
received replacement sub-almanacs, to acquire satellites and take measurements. At 508, 
the MS calculates a coarse position using the sub-almanacs. At 510, the MS transmits the 

20 coarse position to the network along with a list of the satellites used in the solution and 
the associated sub-almanac identifications. The network calculates the final position 
solution at 512 from the information transmitted from the MS. 
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In one embodiment, the method for computing and correcting the coarse location 
includes the following equations. 

The MS computes the coarse location using: 



!±x a ={g t a R- x G a Yg t a R-^p, 



Equation 1 

With: 



Ax, = 



{i-r„J-C'B u \ 



G . = 



Equations 2, 3, 4, and 

10 Equations 5, 6, and 7 

Where: 



A Pu =Pu~ Pi Pu = ^ * k ~ r u ]~ c ■ B <a + c ■ B u 



r„ = 



r « = 



c = 



Pi = 

R = 

B u = 
n = 



Reference location vector in ECEF coordinates (3x1) 

Receiver Location vector at reception time in ECEF coordinates (3x1) 

Speed of light in vacuum 

Almanac-predicted i* satellite location at transmit time in ECEF 

Almanac-predicted ith satellite clock bias 
Measured nseudoranee for ith satellite 
Measurement weighting matrix fnxn) 
Receiver clock bias 

Number of satellites used in location solution 
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The BS corrects coarse location using: 



Ax E = Ax A +P-Ap 
Equation 8 

With 

PiE = X m • [n E ~ K ] ~ c • B iE + c ■ B u 
10 Equation 9 

A 

Ap =P iE -Pu 
Equation 10 

M E = It - r u J - c ■ B u ] P = {g t e R~ 1 G e )G T E R- 
Equations 11, 12, 13 and 



G E = 



Til 



15 



1 - r 'E r u 
iE ~ |_ - 

\ r iE-r u 



Equation 14, where: 



r iE = 



B iE = 



Ephemeris-predicted i th satellite position at transmit time (nxl) 
Ephemeris-predicted i* satellite clock bias 



Simulations were run to evaluate the accuracy and the sub-almanac data 
transmission requirements of the coarse location positioning methods described herein. 
The simulations were run using archived sub-almanac and ephemeris data. The 
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simulations were started with a random reference position 10 kilometers from the true 
position of the unit under test and an up-to-date almanac. Time was moved one week per 
iteration, and the initial almanac was used to calculate a coarse position at the MS. The 
coarse position and the sub-almanac identifications for the satellites used in the solution 
5 were made available to a simulated location server. The location server, using the 

ephemeris data appropriate for the simulated week, then calculated the corrected position. 
At each step, if a satellite's sub-almanac projected an out-of-tolerance range error, a sub- 
almanac current to that simulated week replaced that satellite's sub-almanac and the 
replacement was noted. The results of numerous simulations showed that over a 26 week 
10 time period, accuracies averaged between 15 and 25 meters. Over the 26 weeks in the 
simulation, between 6 and 12 satellites required a single sub-almanac replacement. 

Results of two simulations are presented in the following figures. The results of 
one simulation are shown in Figures 6A-6H, which illustrate results of a simulation 
using a 20 kilometer range error replacement threshold. 

15 Figure 6A illustrates calculated horizontal position error for the simulation. 

Figure 6B illustrates the number of sub-almanacs that required replacement in the 
simulation. Figure 6C shows the range error estimate for the satellite with the 
pseudorandom number code (PRN) 3. Figure 6D shows the range error estimate for the 
satellite with the pseudorandom number code (PRN) 23. Figure 6E shows the range 

20 error estimate for the satellite with the pseudorandom number code (PRN) 17. Figure 6F 
shows the range error estimate for the satellite with the pseudorandom number code 
(PRN) 26. Figure 6G shows the range error estimate for the satellite with the 
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pseudorandom number code (PRN) 6. Figure 6H shows the range error estimate for the 
satellite with the pseudorandom number code (PRN) 21. 

The results of another simulation are shown in Figures 7A-7H, which illustrate 
results of a simulation using a 50 kilometer range error replacement threshold. 

5 Figure 7A illustrates calculated horizontal position error for the simulation. 

Figure 7B illustrates the number of sub-almanacs that required replacement in the 
simulation. Figure 7C shows the range error estimate for the satellite with the 
pseudorandom number code (PRN) 27. Figure 7D shows the range error estimate for the 
satellite with the pseudorandom number code (PRN) 31 . Figure 7E shows the range 
1 0 error estimate for the satellite with the pseudorandom number code (PRN) 8. Figure 7F 
shows the range error estimate for the satellite with the pseudorandom number code 
(PRN) 23. Figure 7G shows the range error estimate for the satellite with the 
pseudorandom number code (PRN) 29. Figure 7H shows the range error estimate for the 
satellite with the pseudorandom number code (PRN) 7. 

15 To illustrate the downlink bandwidth savings using embodiments of coarse 

positioning as described herein when compared to a typical acquisition assistance 
method, two scenarios with their respective bit usages are provided below. Uplink 
messaging is comparable between the two techniques. 

One scenario is turn-by-turn directions. The list below shows particulars used in 
20 the comparison between conventional acquisition assistance and coarse location 

positioning as described herein. Table 1 shows that a reduced number of assistance bits 
are required when using embodiments of coarse location positioning as described herein. 
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Number of users: 10 

Average length of the session; 20 minutes 

Usable age for acquisition assistance data: 4 minutes 

For the coarse location positioning, assume each satellite gets a new sub-almanac during 
5 the session (very pessimistic; this would only occur if the user had not performed an LCS 
session during the previous 3 to 20 weeks) 

The comparable downlink bandwidth used is: 

10 TABLE 1 



Number of Satellites 
in the Solution 


Total Assistance Bits 
Using Acquisition Assistance 


Total Assistance Bits 
Using Coarse Location 
Positioning 


4 


12,550 


7,660 


6 


18,150 


11,420 


8 


23,750 


15,180 



Another scenario is asset tracking in which the assets are taxi cabs. The list below 

shows particulars used in the comparison between conventional acquisition assistance 

and coarse location positioning as described herein. Table 2 shows that a reduced 

15 number of assistance bits are required when using embodiments of coarse location 

positioning as described herein. 

Number of taxis: 15 
Average length of the shift: 8 hours 
Usable age for acquisition assistance data: 4 minutes 
20 For coarse location positioning, assume each satellite gets a new sub-almanac during the 
shift, (again, very pessimistic; this would only occur if the user had not performed an 
LCS session during the previous 3 to 20 weeks) 

The comparable downlink bandwidth used in this scenario is: 

25 

TABLE 2 



Number of Satellites 
in the solution 


Total Assistance Bits 
Using Acquisition Assistance 


Total Assistance Bits 
Using Coarse Location 
Positioning 


4 


451,800 


1 1 ,490 


6 


653,400 


17,130 


8 


855,000 


22,770 
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Various embodiments of the invention have been described with reference to the 
figures, equations, and tables, none of which are intended to be limiting. Within the 
scope of the invention, as claimed below, are multiple alternatives not specifically 
5 described. For example, the order of operations in the illustrated methods may be 

changed without departing from the scope of the invention. In addition, the components 
and their respective functions as shown may be rearranged or reassigned without 
departing from the scope of the invention. 
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